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The scaling behavior of factorial moments of the differences in multiplicities between 
well separated bins in heavy-ion collisions is proposed as a probe of quark-hadron phase 
transition. The method takes into account some of the physical features of nuclear colli- 
sions that cause some difficulty in the application of the usual method. It is shown in the 
Ginzburg-Landau theory that a numerical value 7 of the scaling exponent can be deter- 
mined independent of the parameters in the problem. The universality of 7 characterizes 
quark-hadron phase transition, and can be tested directly by appropriately analyzed data. 



I. INTRODUCTION 



The study of multiplicity fluctuations as a phenomenological manifestation of quark-hadron phase tran- 
sition (PT) has been pursued in recent years in the framework of the Ginzburg-Landau (GL) formalism 
for both the second-order 0- § and first-order @- § PT. In both cases scaling behaviors of the factorial 
moments have been found and are characterized by scaling exponents v. The value of v for second order 
PT is independent of the details of the GL parameters. It therefore provides a distinctive signature for the 
existence of quark-gluon plasma, if its transition to hadrons is of the second order [Q . 

The experimental verification of v has not been carried out so far in heavy-ion collisions, although it 
has been checked to a high degree of accuracy in quantum optics ||. There are a number of reasons for 
the difficulties, which will be described below. The aim of this paper is to devise a method to circumvent 
the obstacles that stand in the way of extracting the signal from the experimental data. In so doing 
we also broaden the scope of the analysis to include aspects of wavelets and correlations, in addition to 
incorporating some evolutionary properties of heavy-ion collisions that are particularly relevant to hadron 
production during PT. 

Let us now examine the difficulties in analyzing multiplicity fluctuations in heavy-ion collisions. The 
factorial moments that have been suggested to quantify the fluctuations are defined by ]9| 

_ (n(n-l)---(n-q + l)) 

where the averages are performed over a distribution, P„, of the multiplicities n in a bin of size 5. Note that 
n must be > q in order for an event to contribute to F q , for which q is usually an integer ranging from 2 to 
5. As 5 is decreased, the average multiplicity (n) in a bin decreases, and may become < g in a hadronic 
collision. Thus F q measures the tail end of the distribution P„, where n 3> (n), i.e. large multiplicity 
fluctuations. Intermittency refers to the power-law dependence of F q on S |l 

F q oc S~ v * , (2) 

a behavior that has been found to be ubiquitous in hadronic and leptonic processes Jlfjfl . However, in nuclear 
collisions the situation is very different. 



In the first place the average multiplicity in an event in nuclear collisions is so high that even in bins of 
small 6 the values of n are large compared to the order q in (jl|) that has been examined experimentally. 
Thus unlike hadronic and leptonic processes the existent nuclear data have not been analyzed to study 
large multiplicity fluctuations, for which q must be increased to values 3> (n) in small bins. The moments 
F q for q = 3,4, and 5 that have been determined are dominated by contributions from the lower-order 
correlations. This point has been emphasized by Sarcevic and collaborators pj]Jl2| , where the nuclear data 
on the cumulants K q are shown to be consistent with zero for q > 3. Whether K q acquire nonvanishing 
values at smaller 8 (so that (n) <C (?) is not known. Thus until the experiments can be improved to render 
the analysis at very small S feasible, a new method must be devised to circumvent this difficulty of extracting 
dynamical information at medium values of q. Our strategy is to consider the distribution of the differences 
in multiplicities between bins and to examine the scaling behavior of its moments. 

Another difficulty associated with the problem of linking multiplicity fluctuations to the dynamics of 
phase transition is that the particles are integrated over the entire duration when quarks turn into hadrons, 
while the system undergoes an expansion. In a second-order PT the fluctuations in hadronization can be 
large relative to the average, but that average is over a short period near T c where hadronization is not 
robust. When integrated over the whole history of the nuclear collision process, such fluctuations may well 
be averaged out, leaving no discernible effect at the detector, which collects all the particles produced in an 
event. This problem is present even if there is no thermalization of the hadrons in the final state, which 
we shall assume in order to focus our investigation here. Our present method is to apply the GL theory in 
increments of time when hadrons are produced near T c , and to integrate the production process over the 
entire duration to yield the measurable multiplicity fluctuations. Our aim is to show that with appropriate 
care in treating the moments certain scaling behavior persists and characterizes the dynamics of phase 
transition. 

Our result reveals a new scaling exponent 7, different from the one, v, found in Refs. @,|. It is not a 
revision of v but a new exponent, since different quantities are investigated. Independent of the theoretical 
considerations underlying this work, the proposed moments can be determined experimentally. Nuclear data 
should be analyzed in the way suggested, even if PT is not an issue. If the scaling behavior is found, but 
the scaling exponent does not agree with the predicted 7, it would not only imply that there has not been 
a PT of the GL type, but also present a numerical objective for a successful hadronization model to attain 
for heavy-ion collisions. 



II. MULTIPLICITY DIFFERENCE CORRELATORS 



To overcome the problem of high multiplicity per bin in heavy-ion collisions we introduce the factorial 
moments that we shall refer to as multiplicity difference correlators (MDC). They are a form of hybrid of 
factorial correlators Q and wavelets |l3| . MDC are not as elaborate as either one of the two separately, but 
are simpler combinations of the two, possessing the virtues of both. 

Let us start with a brief look at the application of wavelet analysis to multiparticle physics [ jLUJli and 



astrophysics [ p_5| . Given any event, the structure of a spatial pattern, exemplified here by a one-dimensional 
function g(x), can be analyzed by a multiresolution decomposition, using the basis functions ip^ k (x) such 
that the wavelet coefficients are 

W jk = t dx i>f k {x) g(x) , (3) 



where L is the domain of g(x) being analyzed. This is like a Fourier transform, except that instead of the 
exponential factor exp(inirx/ L), one uses the Haar wavelets 

^f k {x)=^ H {Tx-k) , (4) 

where 

1, < x < 1/2 

ij) H {x) = { -1, 1/2 < x < 1 (5) 
0, elsewhere . 
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The basis function tpK(x) has a scale index j and a translation index k, which enable the wavelet analysis 
to identify the location and scale of an arbitrary fluctuation of g (x) in terms of Wjk- 

While (Q) and other transforms like it are very powerful and have various virtues, such as their invertibility, 
they offer more than what is needed at this point in our search for a measure of the multiplicity fluctuations 
in heavy-ion collisions that can convey the signature of phase transition. A complete set of Wjk records all 
the information in the spatial pattern in an event, whereas we look for an efficient way of capturing some 
general features averaged over all events. 

The Haar wavelets, however, contain some features that we regard as important ingredients for an im- 
provement of the usual factorial moments and correlators |p|, ^0|Jl7| . For a fixed set of the indices j and fc, 
we have 



^f h (x) =+ 1, for xex+ = i^x\ kjV <x<(k+~^ /V 



I. (ur.r ; z = \.v\ (fc+i) rl> ' , • </,• -! 1.1/2'' } (()) 



Wjk=n+-n-, n±= / dxg(x) . (7) 

J x± 

It is the difference between the contributions in the neighboring bins, x+ and X-. In some wavelet analyses 
moments that involve sums over k are studied to reveal scaling behaviors in j |L5|]. That can be done 
for one event, sometimes referred to as horizontal analysis. For our purpose in what follows, we prefer 
to emphasize first the vertical analysis; i. e., we average over all events for fixed bins. For multiparticle 
production g(x) would be the particle density and n± the multiplicities in the bins x±. Studying the 
difference n + — n_ of neighboring multiplicities is a way to overcome the problem of high multiplicity per 
bin without abandoning the focus on multiplicity fluctuations. 

Once we consider multiplicity differences there is no reason why we should restrict the two bins to only 
the neighboring ones. Let A be the distance between two bins, each of size 5, and let n\ and n 2 be the 
multiplicities in those bins in a given event. Define m to be the multiplicity difference, m = n\ — n 2 . We 
shall be interested in the distribution Q m (A,S) after sampling over many events at fixed A and S. The 
factorial moments we shall study are the MDC 

T q = f g = ^2m(m-l)---(m-q+l)Q m . (8) 

^ 1 m=q 

They are not the same as the Bialas-Peschanski correlators, which are normalized products of factorial 
moments of multiplicities in two bins, averaged over all events, 

p _ (ni (ni - 1) ■ ■ ■ (ni - qi + 1) n 2 {n 2 - !)■■■ (n 2 - g2 + 1)) (Q , 

9192 / \9l / \92 ' V*) 

but the two are similar to the extent that n\ and n 2 are the multiplicities in the two bins separated by A. 
F qiq2 have been found to depend on A, but not on S [0, but T q as defined in (||) will depend on both. T q 
differ from F q defined in ([[]) in that they are the moments of the multiplicity difference distribution (MDD) 
Q m (A,<5), which is a generalization of the usual multiplicity distribution in a way that incorporates the 
virtues of both wavelets and correlators. 

Before we enter into the theoretical description of the MDD Q m in the following sections, we note that 
the experimental determination of T q for hadronic and nuclear collisions should be performed independent 
of theory. Their dependences on A and 8 will pose a challenge to any model of such collisions. Our concern 
in this paper is to make a theoretical prediction of what should be observed when there is a quark-hadron 
PT. But if there is no PT, the properties of T q will remain as valuable features of multiparticle production 
that a good model of soft interaction must explain. 
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III. STATISTICAL AND DYNAMICAL FLUCTUATIONS 



Although our aim is to study the nature of the fluctuations due to quark-hadron PT in heavy-ion colli- 
sions, we begin with a formulation that is more generally valid for any hadronic or nuclear collisions. The 
multiplicities in the two bins discussed in the preceding section fluctuate both statistically and dynamically. 
Focusing on just the statistical part first, and using, as usual, the Poisson distribution V ni for it, we have 
the statistical MDD 

(10) 

ni ,n 2 

where si and s 2 are the average multiplicities in the two bins, and 

-PnM) = —^re- s > . (11) 

Note that, unlike the usual multiplicity, m can be both positive and negative. In fact, the sums in ( |Io| ) can 
be analytically performed, yielding 

Pn (S1,S 2 ) = ( Sl / S2 y n/2 Irn (2y/^) e^ 1 ^ 2 , (12) 

where I rn is the modified Bessel function. 

In the following we shall not assign any intrinsic properties to the two bins and consider only the absolute 
difference between their multiplicities, i. e. 

m = |m — n 2 \ ■ (13) 
Since I m is symmetric under m <-*• — m, we have from ( |l2|) with m > 

P m (si,s 2 ) = cosh( _ln-i W m (2V^)e- Sl ^ 2 (2-<$mo) . (14) 

There is a discontinuity at m = because for all m > there is a reflection of — m in (|l2|) to +m. P m is 
properly normalized to J2m=o Prn = 1 ■ 

If s\ and s 2 are large, then V ni {si) may be well approximated by Gaussian distributions, and P rn would 
also be Gaussian with a width proportional to (sis-2) 1 ^ 4 . It is because of this reduced width that we consider 
MDD: as discussed in Sec. 1, we need smaller values of m to render lower-order moments effective in 
measuring the fluctuations. If s\ and s 2 are small, then P m becomes Poissonian also. For that reason we 
shall consider factorial moments of P m , since the statistical fluctuations can thereby be filtered out |J. In 
both experimental analysis and theoretical consideration the bin width 5 is to be varied so that s, will range 
over both large and small values; hence, no approximation of P m will be made. 

We now introduce the dynamical component of the fluctuations. Denoting it by D (s±, s 2 , A, 5) we have 
for the observable MDD 



Q m (A, S) = J dsids 2 P m (si,s 2 ) D (si, s 2 , 



A, 5) . (15) 



In essence this is a double Poisson transform of the dynamical D, which is a generalization of the formalism 
for photon counting in quantum optics, later adapted by BP for particle production ||. It is this distribution 
that we have proposed to study by use of the MDC T q , defined in m). Our aim is to examine the scaling 
properties of T q and extract universal features that are characteristic of the dynamics of the problem. 

There are many directions in which one can pursue from here. For the dynamical distribution D one 
can consider the mathematical a model [fill, or the more physical models such as the Fritiof model Ji"9| |, 
VENUS [^0| and ECCO |2l[]. Alternatively, one may want to emphasize the large-scale space-time structure 
by considering the small A behavior in an interferometry type of analysis p2fl . For us in this paper we want 
to consider the opposite, namely: the large A behavior where the global size of the particle emitting volume 
is unimportant. The usual short-range correlation in rapidity in low-px multiparticle production JlOj] would 
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also be not important, if A is sufficiently large, but not large enough to cause kinematical constraint. Our 
purpose is to go to a region where fluctuations due to the dynamics of PT are the only ones that need to be 
taken into account. 

More specifically, we shall identify D (si, S2, A, 5) with the Boltzmann factor, exp(— F), in the Ginzburg- 
Landau theory, in which si and S2 will be functions of the order parameter. In fact, for two identical bins 
at large A apart in two regions of the expanding system that have similar spatial and temporal properties, 
we may set s\ = s 2 , ignore A and rewrite ( |l5| ) as 

Q m {5) = J dsP m (s)D(s,S) . (16) 

There remains a complication arising from the property that PT occurs over an extended period of time and 
that the detected hadrons are the result of an integration over that period. That is the subject we turn to 
in the next section. 



IV. HADRONIZATION IN GINZBURG-LANDAU THEORY 



The conventional view of the physical system in a heavy-ion collision at very high energy is that a cylinder 
of locally thermalized partons expands as a fluid, mainly in the longitudinal direction, but also in the radial 
direction at a slower rate. If the colliding nuclei are massive enough, and the incident energy high enough, 
the temperature in the interior of the cylinder may be higher than the critical temperature T c for quark- 
hadron PT, which we shall assume to be second order. Due to the transverse expansion the temperature 
T decreases with increasing radius, at least initially and for the most part of the lifetime of the system. 
Thus hadronization takes place mainly on the surface of the cylinder where T T c . Being a second-order 
PT there is no mixed phase where quarks and hadrons coexist. We assume that there is no thermalized 
hadronic phase surrounding the partonic cylinder, so the hadrons formed on the surface move in free flow 
to the detector. With these simplifying assumptions, which are not unrealistic at extremely high collision 
energies, we can then focus on the issue of relating the hadronization process on the surface to the hadron 
multiplicity collected by the detector. 

Consider first just one bin, which occupies 5r] 8(p in pseudorapidity rj and azimuthal angle ip. We define S 2 
to be its area. Such an area selected by an analyst of the data corresponds to a similar area on the surface 
of the cylinder. Let the hadronization time be th, which is of order 1 fm/c: it is the average time for the 
formation of one hadron on the surface. During that time we use the GL description of PT to specify the 
probability D(s, 6) that s hadrons are created in the area S 2 . The GL free energy, being time independent, 
does not track the time evolution of the system. It is also not equipped to describe the spatial variations on 
the surface, nor need it be. As a mean field theory, it is concerned with the probability for PT near T c as a 
function of the order parameter <fi. Let the two-dimensional coordinates on the surface be labeled by z, then 
the GL free energy is p3] 



F[4>] = / dz 
Js 2 



a|<Mz)| 2 + fo|0(z)| 4 + c|cW^| 2 . (17) 



where a, b, and c are GL parameters that depend on T. For hadronization within the bin of size 5 2 it is 
only necessary to integrate z over that area. For T ~ T c the GL theory requires that a < 0, and b > 0. 
We have found in [QJ3) that for small bins the third term in (|l7]) does not have any significant effect on the 
multiplicity fluctuations, so we shall set c = 0, as it has been done in all previous work Jl[- Furthermore, 
we make the approximation that (f>(z) is constant in S 2 so that 



(a\4 



&I0I 4 ) , (18) 



Assuming that this is valid for any local area on the cylindrical surface, the same free energy is to be used 
for both bins of P m in ([l5|), resulting in the same average multiplicity s\ and S2- We now must specify the 
relationship between s = s\ = S2 and |</>| 2 . 
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As discussed extensively in the square of the order parameter |(/>| 2 is the hadron density p, which, in 
the absence of fluctuations around the minimum of (18) at \(j>o\ 2 , is zero for T > T c but positive for T < T c . 
Thus the average multiplicity in a bin is 

n = 6 2 \ ( f ) \ 2 . (19) 

It should be borne in mind, however, that this is the average multiplicity during the hadronization time 
th, when the system is momentarily regarded as static, and the GL consideration is applied to describe 
the formation of hadrons from quarks at T ~ T c , It is not the average multiplicity in S 2 registered by the 
detector, since the experimental measurement integrates the hadronization process over the entire lifetime T 
of the whole parton system, during which partons are continuously converted to hadrons. Thus the measured 
average multiplicity in 6 is 

rT 

dtmf , (20) 

where <j){t) is formally the time-dependent order parameter, which has to be varied in a functional-integral 
description of Q m (S). Indeed, s in (pC|) is the average multiplicity of P m (s) in (|l6|), but no in ( |l9|) is the 
relevant average multiplicity for the GL description during th- Thus to adapt the formalism of the previous 
sections to the PT process in heavy-ion collisions, it is necessary to modify ( |l6| ) into a functional integral 

g m (^r)-Z- i y^P m ( ( 5 2 r|0| 2 )exp[- ( 5 2 ( a |0| 2 + fe|0| 4 )] , (21) 

where Vtfi = 7r<i|0| 2 , Z = j T><f>exp [— F((f>)] and the integral in ( po|) has been discretized into r = T/th 
segments. In each segment </> is spatially and temporally constant in S 2 and is itself integrated over the whole 
complex plane. For collisions of large nuclei r may be a large number. Herein lies the crux of the problem: 
PT at T ~ T c gives sparingly few hadrons within 5 2 in any time interval around t^, but the detected number 
in a bin at the end of the collision process is roughly r times as many. Our aim to find some universal feature 
of the problem that is essentially independent of a, b, 5 and r. 



V. SCALING BEHAVIOR 



We now proceed to determine the multiplicity difference correlator MDC, which in the present case of 
large A has no dependence on A and is just the normalized factorial moments J- q of the MDD Q m defined 
in (||). As we have mentioned in general terms earlier, an important reason for studying Q m instead of the 
usual multiplicity distribution P n in a single bin is the largeness of S 2 t \ 4>\ 2 in (pl|), when r is large. In that 
case the bin multiplicity n is high, so factorial moments of low orders are ineffective in extracting genuine 
correlations. Multiplicity difference m is on the average proportional to y/r, thus enabling T q to be more 
effective at low q. 

Setting s% = S2 = s in (fl4|), with 

s = 6 2 t\0\ 2 , (22) 

we have the distribution P m in ( pl| ) 

P rn ( S ) = (2 - 6rn0) Irn(2s)e- 2S , (23) 

where m > 0. Let us simplify ( |2l| ) by using the variable u 2 = S 2 b |</>| 4 , getting 

/•oo 

Q„ 1 (t,x) = Z- 1 du P m (rxu) e xu ~ u2 , (24) 
Jo 

where 
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Z= due xu ~ u , x = \a\5/Vb . (25) 
Jo 

For notational economy r has been redefined as 

7 (26) 



th | a 



It is therefore the r defined earlier in units of \a\. It makes physical sense because \a\ is a measure of how 
readily hadronization can take place in the GL theory. The minimum of F[(j>] in ( |l8| ) is at \4>o\ 2 — | a. | /2b, 
for a < 0, so there would be virtually no hadron condensates apart from fluctuations, if \a\ — * 0, resulting 
in t — > oo. Thus the expanding parton system must drive the surface temperature to below T c , making 
a sufficiently negative and \4>o\ 2 large enough to produce hadrons at a rate just such as to carry away the 
necessary energy to maintain the hydrodynamical flow with T < T c at the surface. The beauty of the GL 
approach is that all the complications of the hydrodynamics of the problem are hidden in a few parameters, 
which would be in the final answer [like r and x in (p4|)] unless we can find observables that are insensitive 
to them. 

To find such observables we now calculate the normalized factorial moments T q of Q m using (^) . We first 
fix r and examine T q as a function of x, which is proportional to i5. In Fig. 1 we show how liLF q depends 
on log 2 x for t = 10 and 2 < q < 6. No simple behavior can be ascribed to the rising and falling of the 
curves. It should be noticed that the increase of T q with increasing S is opposite to the usual behavior 
of the single-bin factorial moments F q which increase with decreasing 5. Thus T q do not have the 

intermittency behavior of BP || . 

The similarity of the dependences of T q on x for the various q values shown in Fig. 1 suggest that we 
should examine, as in ffl, the dependences of T q on T% for 3 < q < 6. That is shown in Fig. 2, where the 
dots are the values of T q for log 2 x — i/2 with i being integers in the range — 12 < i < 4. The straightlincs 
are fits of the linear portions. Clearly, J- q exhibit the power-law behavior 

T q ^T^" , (27) 

which we shall call ^-scaling. It is a scaling behavior that is independent of S and b. The dependence on a 
is so far unknown, since the calculation is done at a fixed r. 

From the slopes of the straight lines in Fig. 2 we show in Fig. 3 the dependence of (3 q on q. It can be well 
fitted by the formula 

P q = {q-lV , (28) 

where 7 = 1.099. We use the symbol 7 here for the scaling exponent, instead of v, which we have used 
previously for a similarly defined quantity [as in (^8|)] for F q . For comparison, we recall that [|] 

v = 1.304 . (29) 

Thus the scaling exponent 7 for MDC is significantly lower. Lower value of that exponent means larger 
fluctuations. 

So far the result is for r = 10 only. To see the dependences on r, we have repeated the calculation for a 
range of r values. Scaling behavior as in ( p7j ) has been found in each case, and ( p8| ) is also well satisfied. 
Fig. 4 shows the dependence of 7 on r. Evidently, it is nearly constant for 3 < r < 30 with the value 

7 = 1.09 ±0.02 . (30) 

Thus, the result of this study is embodied in just one number, 7. It is independent of a, 6, 5, and r, provided 
that a is negative to allow hadronization. 

The universality of 7 is remarkable and should be checked experimentally. If a signature of quark-hadron 
PT depends on the details of the heavy-ion collisions, such as nuclear sizes, collision energy, transverse energy, 
etc., even after they have passed the thresholds for the creation of quark-gluon plasma, such a signature is 
likely to be sensitive to the theoretical model used. Here 7 is independent of such details; it depends only 
on the validity of the GL description of PT for the present problem. The MDD Q m can readily be measured 
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experimentally, and the moments T q directly determined as functions of bin size. If the scaling behavior 



(|27|), supplemented by (28), is satisfied with 7 ~ 1.1, then we may interpret the system as having undergone 
a second-order PT describable by the GL theory. 

In the absence of PT we may regard the two bins separated by a large A to be totally uncorrelated. In that 
case Qm(S) can be identified with P m [s(<5)], where, if we vary s in the range 1 < s < 100, the corresponding 
T q without PT can be calculated directly. Again, F-scaling is found, satisfying ( p7| ) and (|2^), but this time 
with 

7 = 1.33 ±0.02 (noPT) . (31) 



This value is sufficiently separated from that of (|30|) derived for PT so that phenomenological distinquisha- 
bility between the two cases should be quite feasible. 



VI. CONCLUSION 



In this work we have solved a number of problems that have obstructed the study of multiplicity fluctuations 
in heavy-ion collisions as a means of finding signatures of quark-hadron phase transition. One problem is 
the large multiplicities even in small bins, for which the usual factorial moments F q fail to reveal distinctive 
features for 3 < q < 6, since events with large fluctuations are submerged by generic events. That is mainly 
an experimental problem where the analysis of the data cannot be pushed to the regions (n) s < q. Another 
problem of more theoretical nature is that the application of the Ginzburg-Landau theory of PT needs 
special tailoring for a system whose lifetime is finite, but long compared to the transition time for individual 
hadrons, and whose observables are integrated over that time. We have overcome both of these problems by 
showing the effectiveness of studying the fluctuations of multiplicity differences. 

We have started with wavelet analysis and found that it generates more information than can easily be 
filtered to yield a succinct signature of PT. However, we extracted a simple feature of the wavelets and 
considered the MDD, Q m , involving two bins separated by a distance A. Although A can be any value, we 
have considered only large A in order to apply the simplest description of PT by the GL theory. The time 
integration problem of the detectable multiplicities is handled at the expense of an extra parameter t, which 
cannot be specified without a hydrodynamical model lying outside the scope of this treatment. The goal has 
then been to find a measure of PT that is independent of the unknown parameters in the problem. 

That goals was achieved by the discovery that the factorial moments T q of Q m satisfy a scaling behavior 
that is characterized by a number 7 ~ 1.1. It is independent of the details of the dynamics, except that a 
PT occurs in a way describable by the the GL theory. Thus 7 is a universal constant for the problem. We 
call 7 a scaling exponent, but it has absolutely no connection with the critical exponents of the conventional 
critical phenomena. There is no need to know the temperature, which is not measurable, or the critical 
T c . The moments T q can directly be determined from the data, and can therefore be checked for F-scaling, 
independent of any theoretical input. Whatever experimental value the exponent 7 turns out to be would 
be of great interest, since a value different from 1.1 would demand an explanation. If it approaches 1.1 as 
energy or nuclear size is increased, it would suggest an approach to the condition required for the formation 
of quark-gluon plasma. 

Experimentally, it should be easy to vary A and check not only how T q themselves depend on A, but 
also whether and how the scaling behavior is affected. The A dependence has not been investigated here. 
Theoretically, there are many other challenges that also await undertaking. A verification of 7 = 1.1 will 
undoubtedly stimulate an upsurge of interest in the study of multiplicity fluctuations in particle and nuclear 
collisions. 
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FIGURE CAPTIONS 



Fig. 1 Factorial moments T q versus x in log-log plot for various values of q at r = 10. 

Fig. 2 Factorial moments T q versus Ti in log-log plot for various values of q at r = 10. Dots are calculated 
values at x = 2 4 / 2 for —12 < i < 4. Straight lines are fits of the linear portions of the dots. 

Fig. 3 Dots are the slopes of the straight lines in Fig. 2, plotted against q— 1 in log-log plot. The straight 
line is the best fit of the dots, whose slope is 7, defined in Eq. (28). 

Fig. 4 Solid line shows the dependence of 7 on r; dotted line is at 7 = 1.09 to guide the eye. 
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